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p , 1 Abstract 

,-S^ We consider the supersymmetric extension of the 3-3-1 model with right-handed neutrinos. We 



study the mass spectra in the scalar and pseudoscalar sectors, and for a given set of the input 



parameters, we find that the lightest scalar in the model has a mass of 130 GeV and the lightest 
pseudoscalar has mass of 5 GeV. However, this pseudoscalar decouples from the Z° at high energy 
scales since it is almost a singlet under SU(2)l <S> U(l)y- 

PACS numbers: PACS number(s): 12.60. Jv; 12.60.-i; 14.60.St 



* E-mail address: montero@ift.unesp.br 

t E-mail address: vicente@ift.unesp.br 

■t E-mail address: mcrodriguez@fisica.furg.br 



I. INTRODUCTION 



Models with SU(3)c <8> S77(3)l <S> U(1)n gauge symmetry 



called 3-3-1 models for short) 



QflB 



At low energies they 



are interesting possibilities for the physics at the TeV scale 
coincide with the standard model and some of them give at least partial explanation to some 
fundamental questions that are accommodated but not explained by the standard model. 
For instance, in order to cancel the triangle anomalies, together with asymptotic freedom in 
QCD, the model predicts that the number of generations must be three and only three; (ii) 
the model of Ref. (l| predicts that (g' / g) 2 = sin 2 0^/(1 — 4 sin 2 8 W ), thus there is a Landau 
pole at the energy scale \i at which sm 2 6\y(n) — 1/4. According to recent calculations 



En it 

, |5[ ; (Hi) the quantization of the electric charge |6j 

racter of the electromagnetic interactions [7|] do not depend on the nature of the neutrinos 

i.e., if they are Dirac or Majorana particles; (iv) as a consequence of item ii) above, the 

model possesses perturbative Af = 1 supersymmetry naturally at the /i scale jala|; (v) the 

Peccei-Quinn |ll| symmetry occurs naturally in these models |ljj; (vi) since one generation 

of quarks is treated differently from the others this may be lead to a natural explanation for 

the large mass of the top quarks 12 1. Moreover, if right-handed neutrinos are considered 

transforming non-trivially |3j, 3-3-1 models [l|, |2j can be embedded in a model with 3-4-1 

gauge s ym metry in which leptons transform as (yi,l,u£, I c )l ~ (1,4,0) under each gauge 



factors 



ymn 



The SU(3)l symmetry is possibly the largest symmetry involving the known 
leptons (and 577(4)/, if right-handed neutrinos do really exist). This make 3-3-1 or 3-4-1 
models interesting by their own. 

Models with SU(3) (or 577(4)) electroweak symmetry may have doubly charged vector 
bosons. These sort of bileptons may be detected by measuring several left-right asymmetries 
in M0ller scattering J.14J or, in future lepton-lepton accelerators. It is interesting that one 
these parity non- conserving asymmetry was observed in M0ller scattering for the first time 
only recently and its value is, at present, in agreement with the standard model predic- 
tion [16]. However, in the future those asymmetries in e~e~ colliders (jjj may also be used 
for searching doubly charged particles and a heavy neutral Z'° vector bosons, which is also 
a prediction of those models, may be discovered in e~/i~ colliders |l7|. Singly and doubly 



charged vector bileptons may also be produced in e 7 [IS] or 77 [l£ 



or hadron 120] colliders. 



New heavy quarks are also part of the electroweak quark multiplets in the minimal model 



representation. They are singlets under the standard model SU(2)l ® U(l)y group symme- 
try and in some versions of the 3-3-1 models the electric charges of these heavy quarks are 
different from the usual ones, so that it can be used to distinguish such models from their 
viable competitors. In fact, the pp production of these exotic quarks at the energies of the 
Tevatron have been studied in Ref. [2j| where a lower bound of 250 GeV on their masses 
was found. This sort of models are also predictive with respect to neutrino masses 12211: they 
can implement the large mixing angle MSW solution to the solar neutrino issue |23|, and 
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also the almost bi-maximal mixing matrix in the lepton sector 

Gauge models based on the 3-3-1 gauge symmetry can have several representation con- 
tents |l|, |3J] depending on the embedding of the charge operator in the SU(3)l generators, 

9. = hx 3 -t}\ 8 ) + N I, (1) 

where the °d parameter defines two different representation contents, N denotes the U(1)n 
charge and A3, As are the diagonal generators of SU(3). The supersymmetric version of the 
model of Ref. |l| , with 1? = v3, has already been considered in Refs. (a |9| . In this work we 
build the supersymmetric model based on the representation content of the 3-3-1 model of 
Refs. 3[ which corresponds the the case when d = l/\/3. 

The outline of this paper is as follows. In Sec. |H] we review the non-supersymmetric 
3-3-1 model with right-handed neutrinos, introduce the respective superpartners and the 
superfields. The Lagrangian of this 3-3-1 supersymmetric model, including the soft term, is 
considered in Sec. IHIj the scalar potential and the scalar mass spectra are given in Sec. IIV1 
Finally, the last section is devoted to our conclusions. 

II. THE MODEL 



In this section (Sec. Ill All we review the non-supersymmetric 3-3-1 model with right- 
handed neutrinos of Refs. J3j and add the superpartners (Sec. IIIBjI . The superfields are 
introduced in Sec. Ill CI 



A. The representation content 



Let us first summarize the non-supersymmetric model with the charge operator defined 
by Eq. © with $ = l/y/% |^, i.e.,. 



which implies leptons transforming under the 3-3-1 factors as 



(2) 



J aL 



V 

L 



;i,3,-i/3), 



(3) 



\KJ 



with a = e, h,t and v c a = Cv a , plus the singlets 



&~ (1,1,1) 



(4) 



In the quark sector we have the first two families transforming as antitriplets of SU(3)l 



( d \ 



Q< 



with the respective singlets 



U r 



(3,3*,0), a = 1,2; 



(5) 



V </ 



< L ~ (3*, 1, -2/3) , d c aL , Cl ~ (3*, 1, 1/3) . 



(6) 



The third family transforms as triplet under SU(3)l in such a way that, the triangle 
anomaly cancels out only among the three families and taken into account also the color 
charges. 



Q 



31 



( U \ 

Us 

d 3 
yu'j 



(3,3,1/3), 



(7) 



and their respective singlets 



u c 3L , u't ~ (3*, 1, -2/3), d c 3L ~ (3*, 1, 1/3). 



(8) 



In the scalar sector only two triplets rj ~ (1, 3, —1/3) and p ~ (1, 3, 2/3) are necessary to 
break appropriately the gauge symmetry and also to give the correct mass to all the fermions 



in the model. However, to eliminate flavor changing neutral currents we add an extra scalar 
triplet transforming like 77. 



/ T7° \ 

71 




f x ^ 


7f 


, x = 


X' 


U) 




•x/0 
\X2 ) 



;i,3,-i/3), 



<^ 



;i,3,2/3), 



p --- p 

\pi J 

and we will denote the vacuum expectation values which are different from zero as v 
(r ?1 °)/ v / 2, w = ( X ° 2 )/V2 and u = (p°)/V2. 



(9) 



B. Supersymmetric partners 

Here we will follow the usual notation writing for a given fermion /, the respective 
sfermions by / i.e., I and q denote sleptons and squarks respectively. Then, we have the 
following additional representations 



Q 



o.L 



J aL 



( d \ 

1*0 

U c 

\d' aJ 
(v \ 

"a 
la 



(3,3*,0), Q 3L 



;i,3,-i/3), 



u 3 \ 



\ u ) 



(3,3,1/3), 



(10) 



£l~(i,m) 



U iLi U L 



(3*, 1,-2/3), 4,C L ~(3*,l,l/3), 



(11) 



with a = e, h,t; i = 1,2,3; and a = 1,2. However, when considering quark (or squark) 
singlets of a given charge we will use the notation u c iLl d\ L (tin, d\ L with i(j) = 1, 2, 3. 



The supersymmetric partner of the scalar Higgs fields, the higgsinos, are 



71 
fj~ 



X 



( X°i^ 

x~ 

~o 

\X 2 J 



;i,3, 



■1/3), 



;i,3,2/3), 



(12) 



and the respective extra higgsinos, needed to cancel the chiral anomaly of the higgsinos in 
Eq. (O, are 



if 



P 



( fi'O \ 
71 



x' 



( fl \ 

x' + 

\X2 ) 



;i,3*,i/3), 



~ (1,3*, -2/3), (13) 

\P2 J 

and the corresponding scalar partners denoted by i]',x', p' ', with the same charge assignment 
as in Eq. (HHJ), and with the following VEVs: v' = (rff)/y/2, w' = (x' 2 °)/\^ and u' = 
(p'°)/a/2. This complete the representation content of this supersymmetric model. 

Concerning the gauge bosons and their superpartners, if we denote the gluons by g b the 
respective superparticles, the gluinos, are denoted by A^, with b = 1,...,8; and in the 
electroweak sector we have V b , the gauge boson of SU(3)l, and their gauginos partners A^; 
finally we have the gauge boson of U(1)n, denoted by V, and its supersymmetric partner 
Xb- This is the total number of fields in the minimal supersymmetric extension of the 3-3-1 
model of Refs. [3J. 

C. Superfields 



The superfields formalism is useful in writin g th e Lagrangian which is manifestly invari 



ant under the supersymmetric transformations 



2jj | with fermions and scalars put in chiral 



superfields while the gauge bosons in vector superfields. As usual the superfield of a field 



H- 



will be denoted by |26j]. The chiral superfield of a multiplet is denoted by 

<j> = 4>(x, 9, 9) = 4>(x) + i 9a m 9 d m (p(x) + -99 99 D0(x) 

+ V2 9<S>{x) + -^99 9a m d m <j){x) 
v2 

+ e6F (j> (x), (14) 

while the vector superfield is given by 

V(x,6,9) = -9a m 9V m (x)+t999V(x)-i999V(x) 

+ -9999D{x). (15) 

The fields F and D are auxiliary fields which are needed to close the supersymmetric algebra 
and eventually will be eliminated using their motion equations. 
For fermion superfields we use the notation 

LaLj VaLi QaL, Q3L, U° iL , d° iL , u[ : d£ L . (16) 

For scalar superfields we write: fj, x, p and similar expressions for fj', x', p' an d we must 
change (field) by (field) . 

The vector superfield for the gauge bosons of each factor SU(3)c, SU(3)l and U(1)n 
are denoted by Vc, V ' c\ V, V; and V, respectively, where we have defined Vc = T b V^, 
V = T b V b ; V c = f b V b c , V = f b V b ; T b = \ b /2, f b = -\* b /2 are the generators of triplet 
and antitriplets representations, respectively, and X b are the Gell-Mann matrices. 

III. THE LAGRANGIAN 

The Lagrangian of the model has the following form 

^3315 = &SUSY + £ so ft ' (^) 

where Csusy is the supersymmetric part and £ so ft the soft terms breaking explicitly the 
sup er sy mmet ry. 



A. The supersymmetric Lagrangian 



The supersymmetric part of the Lagrangian is decomposed in the lepton, quark, gauge, 
and the scalar sectors as follow: 



C-susy - ^Lepton + £ Quark + £ Gauge + ^Scalar' 



where 



C 



Leptons 



d 4 e 



f 2gV-^V'f , i g'V'f 



aL^ h aL 



(18) 



(19) 



in the lepton sector, we have omitted the sum over the three lepton family for simplicity, 
and 



^Quarks _ d ® Q»L e 



MgsVc+gV)\ 



QaL + V3L e 



[2(gsVc+gV)+^V'} 



Q 



3t 



+ %^> p °-^\? L + d° iL eto>Vc+Wd! L 



a 



^c^Vc-H-V^c + J { 2g s V c+ ^V^ 



aL 



(20) 



in the quark sector, and we have denoted g s ,g,g' the gauge coupling constants for the 
SU(3)c, SU(3)l, U(1)n factors, respectively. In the gauge sector we have 

1 



£ 



Gauge 



d 2 9 [W C W C + WW + WW 



i / d 2 9 [W C W C + WW + WW ] 



(21) 



where Wq, W e W are fields that can be written as follows 

W (C = 
W c = 



—DDe- 29 ^c De 29sV c 
8g s C 



Sg C 



W, = -^DDD C V', C = 1,2. 



'C 



(22) 



"Escalar 



Finally, in the scalar sector we have 

f) 

+ [ d 2 9W + I d 2 9W, 



fjelW-^r) + xe [29 ^-^' ] X + pe^ + ^^p 



+ ff e^+iv']ff + x ' etov+t-v^ + p> e [2 g v-?£vy 



(23) 



where W is the superpotential. 



B. Superpotential 



The superpotential of the model is decomposed as follows 



Wo W-i - Wo W-x 
2 3' 23 V ; 

^2(^2) having two chiral superfields and W^W^) three superfields. Explicitly we have that 
the term with two superfields is given by: 

r 

W 2 = ^{jioaLrf + fiiaLx] + IXjfyfi + fi x xx 

a=e 

+ P P Pp', (25) 

where Lrj = Lifj^. The term with three superfields is written explicitly as 

T T T T T 

W z = y^y^ \abL aL pi c bL + ^[A 2a e^aLXP + A 3a eL aL r}p] + J^ JZ ^b^L aL L bL p 

a=e b=e a=e a=e b=e 

3 3 3 

+ Y K ^Q^f)'u C iL + K^QiL-q'u'l + Y K 2iQ3LXU C iL + K,' 2 Q 3L x'u'l + Y K 3iQsLp' ' d\ L 
i=l i=l i=l 

2 2 3 2 2 2 3 

+ Y K 3aQ^Lpd C a L + Y Y ^aiQaLVd'iL + YY ^QaLV^L + YY K ^iQ a Lpu\ L 
o=l o=l i=l o=l (3=1 a=l i=l 

2 2 3 2 2 

+ Y K 'laQ«LPu£ + YY K ^iQ*LX& lL + Y Y K 5apQaLXdp L 

a=l a=l i=l a=l /3=1 

2 2 2 3 2 3 

+ Y Y Y e f^^QaLQ(SLQyL + h^pxri + f^pxv + Y Y Y tW%Ld 

a=l (3=1 7=1 i=l (3=1 j=l 

3 2 3 3 3 3 3 

+ Y Y ^ c lL d'i L u ic L + YYY ^kd c iL d c jL u c kL + YY ^AA 

1=1 (3=1 1=1 J = l J = l 1=1 J = l 

223 22 t23 

+ YYY ^(3id'^ L d'p L U C iL + YY &a(3d'a L dp L u'Z + YYY ^a(xjL aL Q aL d C jL 
a=l (3=1 i=l a=l (3=1 a=e a=l i=l 

r 2 2 

+ /;/;/; is,aapL aL Q aL d'^ L . (26) 

a=e a=l /3=1 

As usual it is necessary to introduce the so called soft terms that break the supersymmetry 
explicitly. 



(3L u jL 



1 jL u 'L 



C. The soft terms 



The soft terms can be written as 



£ soft - C gmt + ^Scalar + CsMT ' ( 27 ) 



where 



6=1 6=1 

+ m'\ B \ B + H.c.}, (28) 

give mass to the boson superpartners and 

^Scalar = _m ^ tr 7 _ m pP f P - m\rfx ~ mffirf - m 2 p ,p' ] p' - m^x'V 

+ [k 1 e ijk p i XjVk + K e ijkp'iXjVk + H - C -}' ( 29 ) 

in order to give mass to the scalars, we have omitting the sum upon repeated indices, 
i,j, k = 1, 2, 3 . Finally, we have to add 

- Csmt = m\l\ L L aL + m 2 J c JJ c aL + m 2 Q3 Q\ L Q 3L + m 2 QaL Q ] aL Q aL + m 2 u .uf L u c iL + m 2 d .df L d c iL 
+ m 2 u ,ul ] u'l + m 2 d ,j!p L d'p L + [M 2 L aL r) n + M' 2 L aL x ] + e la L aL p'l c L + e 2a eL aL xp 
+ e 3a eL aL r]p + e 4ab eL aL L bL p + QuQ%Lriu c iL + Q^Qzl^u'I + Q2iQzLX'u c iL + q' 2 Qzlx'u'1 

+ Q3aiQaLVdiL + QsapQaLV^L + Q4aiQaLpU C iL + Q^QaLPu'l + Q5iQzLp'd C iL 

+ Q^aQsLP d C aL + QeaiQaLXd'iL + Qf> a pQaLXdp L + QIol^QollQ^lQ^L + v li/3jdi L dp L U < j L 

+ v 2i pd c iL d^ L u[ + v 3ijk d c iL d c jL u c kL + v Aij d c iL d c jL u[ + v 5api d^ L d^ L u c iL + v Ga ^ L d!^ L u'l 
+ v 7aa jL aL Q aL d^ L + v% aa pL aL Q aL dp L + iJ.C], (30) 

in order to give appropriate masses to the sfermions. 

IV. THE SCALAR POTENTIAL 

In the present model the scalar potential is written as 

V 331 = V F + V D + V soiu (31) 



10 



where 



V F 



+ 



-c f 


in 












D 

ijk 


pi) , . fi 

y^ + -j^jkPjXk 


2 
+ 


Px 1 , fz 


2 
+ 


Pv 


1/2 / / 


2 
+ 


Px 1/2 / / 


2 
+ 


P P 
2^ 



yft + Y^jkXjVk 



1/2 / / 
Pi + -j e ^XjVk 



(32) 



and 



Vi 



D 



+ 



-£c = lp B D fl + DD) = ^(-rfv + »/V ~ X ] X + X ] X + 2p f p - 2pV) 2 

£faM,-t& - n?KM + x\\\ 3 x 3 - x'lK^x' + Mpj + pUtpW 



(33) 



finally, 



Koft = -Atft = mffiri + m 2 p p j p + m^X + ra^V 
+ rnfypftpf + m 2 x ,x n x' ~ ^jk(hpiXjVk + Ktix'jVk 
+ H.c), 



(34) 



where we have used the scalar multiplets in Eqs. (0) and (JT3J). 

With Eqs. (|32* jl -(j53 )l we can work out the mass spectra of the scalar and pseudoscalar 
fields by making the usual shift X° — > -4= (vx + ^x + z'-Px)- The analysis is similar to that 
of Ref. 8J and we will not write the constraints equation, etc. 

By using as input the following values for the parameters: sin 2 ^ = 0.2314, g = 0.6532, 
g' = 1.1466; f 2 = 2, f 2 = 10~ 3 ; k t = k[ = 10 GeV; \x n = p p = p x = -10 3 GeV; m v = 15 
GeV, m n = 10 GeV. m n = 244.99 GeV; m v „ = ttkl = 10 3 GeV and m n , = 13 GeV, we 



obtain the masses 



) ""X2 



X2 



m 1 « 1702, m 2 « 1449, m 3 rj 387, 
m 4 « 380, m 5 « 361, m 6 « 130, 



(35) 



for the scalar sector (all masses are in GeV). Note that the lightest neutral scalar is heavier 
than the lower limit of the Higgs scalar of the standard model, i.e., m# ~ 114 GeV. For the 
pseudoscalar sector we obtain 



Mi « 1702, M 2 « 1449, M 3 « 363, 
M 4 « 5, M 5 = 0, M 6 = 0, 



(36) 
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only the two massless pseudoscalars are exact values, i.e., there are two Goldstone bosons 
as it should be. Notice that there is a light pseudoscalar. Although the values of 
the masses above are only an exercise, and it is possible that other values of the in- 
put parameters would give another set of masses, we would like to point out the fol- 
lowing. In the basis (F Vl ,F p ,F X2 ,F r} ',F p /,F x >) the lightest pseudoscalar is given by 
(0.0099,0.0012,0.7070,0.0071,0.0170,0.7070), hence we see that it is mainly F X2 and F^. 
So, we need to investigate the couplings of these pseudoscalars with the Z°. 

In the (W3, Ws, B) basis we have the mass square of the real vector bosons given by: 



/ 



g_ 

4 



v 2 + u 2 



\ 



M v2 - U2 ) 


-f (V 2 + 2U 2 ) 


\ (V 2 + U 2 + 4W 2 ) 


-^=( V 2-2U 2 -2W 2 




^£(V 2 + 4U 2 + W 2 ), 



\ 



(37) 



where we have defined V 2 



v 2 + v' 2 , U 2 = u 2 + u' 2 and W 2 = w 2 + w' 2 , and 

2 -2, 



t z 



sin'fl 



IV 



I sin 2 6 



IV 



The eigenstates of Eq. (|3T|) are 



A 



Vs 



4t 2 + 3 



tw* 



Vs 



W 8 + B 



for the photon, and 



and 



3t 



At 2 + 15t 2 + 9 



t 2 + 3 
3t 



W 3 - -^W 8 + B 
v3 




(38) 



(39) 



(40) 



(41) 



for the Z° and Z ', we have neglected the mixing among Z° and Z ', so that M\jMy V w 
(3 + 4t 2 )/(3 + 1 2 ) = 1/ cos 2 9w- With this at hand, we were able to verify that the couplings 
of F X2 , F x > 2 are given by the usual vertex of the Higgs scalar in the standard model times a 
factor proportional to (W/v w )(U/W) A or (W/v w )(V/W) 4 , where v w ph 246 GeV, and these 
couplings go to zero as W — > 00. This behavior is expected since both x°5 X2 are singlets of 
SU{2) L ® U(1)y and do not couple to the Z° in this limit. 

For completeness we show that the lightest scalar, 

in the basis (H Vl , H p , H X2 , H^ , H p , H x /), is written as 

(-0.0581,-0.9775,0.0610,-0.0394,-0.0592,0.1800), that is, it is mainly H p which 
transforms as doublet under SU(2) L ® U(1)y- 
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V. CONCLUSIONS 

We have put forward a supersymmetric version of the 3-3-1 model of Refs. J3j which in- 
cludes right-handed neutrinos transforming non-trivially under SU(3)l. This sort of models 
is an interesting possibility since neutrinos are massive particles and right-handed neutrinos 
should to be added eventually to any extension of the standard model. 

Concerning the scalar and pseudoscalar mass spectra we have found two different situa- 
tions: for the scalar sector we were able to find all the Higgs scalars with masses above the 
Z° mass and above the lower limit of the standard model Higgs boson obtained by LEP: 
rriH ~ 114 GeV; for the pseudoscalars, for the same set of the input parameters, we have 
found a considerably light one (M 4 = 5 GeV) which in principle can bring some problems. 
However, a carefully analysis have shown that the mass eigenstate corresponding to M4 is 
mainly formed by the symmetry eigenstates F X2 and F x > 2 , and studying the couplings of 
these pseudoscalars with the Z° we have noted that they vanish in the limit where w and w' 
go to infinity i.e., both pseudoscalars decouple from Z° in this limit since they are singlets 
under SU(3) L ® U{1) Y . 

In the other supersymmetric 3-3-1 model $], the proton decay modes are p — > K + u ll 
and p — > K + e T fi ± u T but in the present one, only the mode p — > is^v^ is possible. It means 
that the constraints coming from SuperKamiokande data on p — ► K + u e , which implies that 



Tp > 10 years in this decay mode 



30j, are avoided. 

However, there are higher-dimension (nonrenormalizable) operators, that arise from new 
physics at some scale A. For instance, there are dimension-5 operators that violate baryon or 
lepton number, that are allowed by the gauge invariance, and that contribute to the proton 
decay unless they are sufficiently suppressed. In the context of the MSSM we have J3l| 

^fQiQAU + %W*t (42) 

These terms contribute to the proton decay with diagrams at 1-loop level involving gauginos 
and gluinos, known as dressing diagrams. These are the dangerous terms that induce the 
decay mode p — > K + v e . This channel is enough to exclude the minimal supersymmetric 



SU(5) model 



since the SuperKamiokande data 



H- 



In the present model dimension-5 operators that violate lepton and baryon number, that 
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are allowed by the 3-3-1 symmetry, may arise in unification theories J33j, and are given by 



K 1 K 2 k-3 k- 4 

^aaP/f A wA A n ^ajjrfc fe-c-cjc , ^abi ( f t \(f\ jc\ i ai J 



c 



C fj 7 Q 

+ ^(L a r/)(L^) + ^(L a ^ (43) 



while in the SUSY 3-3-1 model of Ref. p these operators are 
i 23 

K aaf3,f A \/A A x , K Q/3 / A ? C wA ;. c n , K abi 



(L a Q a )(Q 3 Q ) + ^(Q Q J C )(Q 3J -) + ^i( La L b )(Q 3 u< 



4 g 

+ -f^4*)(&^) + ^(Lav'KUf,'). (44) 

The suppression of the effective operators in Eqs. (f4*3~j) and (jl^jl in the context of SUSY 3-3-1 
models will be considered elsewhere J35(. 

Finally, we would like to say that the 3-3-1 model with non-supersymmetric right-handed 
neutrinos furnishes a good candidate for self-interacting dark matter (SIDM) since there 
are two Higgs bosons, one scalar, and one pseudoscalar, which have the properties of the 
candidates for dark matter like, stability, neutrality, and that they must not overpopulate 
the universe [34 1 and, in particular for a self-interacting dark matter candidates, they have 
large scattering cross-section and negligible annihilation or dissipation. It means that the 
candidate for SIDM has not to be introduced ad hoc as in other models J36J. This feature 
remains valid in the supersymmetric version that we have developed in this work, but it 
deserves a more careful study. 
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